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Abstract. Using both the Short- and Long-wavelength Spec- 
trometers on board the Infrared Space Observatory (ISO), 
we have obtained infrared spectra of the Herbig Ae systems 
AB Aur and HD 163296. In addition, we obtained ground- 
based N band images of HD 163296. Our results can be sum- 
marized as follows: (1) The main dust components in AB Aur 
are amorphous silicates, iron oxide and PAHs; (2) The circum- 
stellar dust in HD 163296 consists of amorphous silicates, iron 
oxide, water ice and a small fraction of crystalline silicates; (3) 
The infrared fluxes of HD 163296 are dominated by solid state 
features; (4) The colour temperature of the underlying contin- 
uum is much cooler in HD 163296 than in AB Aur, point- 
ing to the existence of a population of very large (mm sized) 
dust grains in HD 163296; (5) The composition and degree 
of crystallization of circumstellar dust are poorly correlated 
with the age of the central star. The processes of crystallization 
and grain growth are also not necessarily coupled. This means 
that either the evolution of circumstellar dust in protoplanetary 
disks happens very rapidly (within a few Myr), or that this evo- 
lution is governed by factors other than stellar mass and age. 

Key words: Circumstellar matter - Emission-line stars - Stars: 
AB Aur, HD 163296 - Pre-main sequence stars - Infrared: 
Stars 



1. Introduction 

Herbig Ae/Be stars are intermediate-mass pre-main sequence 
stars surrounded by disks of gas and dust which might be the 
site of on-going planet formation (see Waters & Waelkens 1998 
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for a recent review). AB Aur (AOVe) and HD 163296 (A3Ve) 
belong to the best studied Herbig stars and are sometimes con- 
sidered prototypical for the entire class. As early as 1933, Mer- 
rill & Burwell remarked upon the similarity of both systems, 
which was confirmed by numerous subsequent authors. Apart 
from the fact that the stellar mass, effective temperature and age 
of AB Aur and HD 163296 are nearly identical (van den An- 
cker et al. 1997, 1998), the similarity between the two systems 
also extends to their circumstellar environment: both AB Aur 
and HD 163296 are examples of relatively isolated star for- 
mation, and are not hindered by confusion with other sources 
(Henning et al. 1998; Di Francesco et al. 1998). 

Both AB Aur and HD 163296 show a rich, variable, 
emission-line spectrum from the ultraviolet to the optical. With 
a few exceptions, most noticeably the observed [O I] emission, 
these lines have been successfully modelled as arising in an 
inhomogeneous stellar wind (Catala et al. 1989; Bohm et al. 
1996; Bouret et al. 1997). Infall of material has been detected 
in both HD 163296 and AB Aur through monitoring of UV and 
optical absorption and emission lines and can be explained by 
the presence of infalling evaporating exocomets (Grady et al. 
1996, 1999). In the infrared, both stars are among the sources 
with the strongest 10 /im silicate feature in emission (Cohen 
1980; Sitko 1981; Sorrell 1990). Sitko et al. (1999) compared 
the 10 fim silicate feature in HD 163296 with that of solar- 
system comets and found a striking resemblance with that of 
comet Hale-Bopp. Neither star shows the 3.29 /im UIR emis- 
sion band present toward many Herbig Ae/Be stars (Brooke et 
al. 1993; Sitko et al. 1999). Basic astrophysical parameters of 
both stars are listed in Table 1 . 

There is strong evidence for the presence of a circumstellar 
disk in both AB Aur and HD 163296. Bjorkman et al. (1995) 
detected a 90° flip of the polarization angle between the opti- 
cal and the ultraviolet in HD 163296, which they interpreted as 
evidence for a flattened, disk-like structure. Mannings & Sar- 
gent (1997) resolved the gaseous disks surrounding AB Aur 
and HD 163296 using CO millimeter wave aperture synthesis 
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imaging. Using continuum measurements at 1.3 mm, the same 
authors also resolved the circumstellar dust disk of HD 163296. 
The AB Aur dust disk has also been resolved in the infrared, 
and shows a surprisingly strong dependence of disk diame- 
ter on wavelength, ranging from 070065 (0.94 AU) at 2.2 /im 
(Millan-Gabet et al. 1999), through 0"24 (35 AU) at 11.7 /im 
to 0"49 (70 AU) at 17.9 /im (Marsh et al. 1995). 

In this paper we present new infrared spectra of AB Aur 
and HD 163296 obtained with the Short- and Long Wavelength 
Spectrometers on board the Infrared Space Observatory (ISO; 
Kessler et al. 1996). We will discuss these spectra and their im- 
plications for the evolution of dust in Herbig systems. In a sub- 
sequent paper (Bouwman et al. 2000), we will describe a model 
for the circumstellar dust disks of AB Aur and HD 163296 and 
apply it to these data. 

2. Observations 

ISO Short Wavelength (2.4-45 /im) Spectrometer (SWS; de 
Graauw et al. 1996) and Long Wavelength (43-197 /im) Spec- 
trometer (LWS; Clegg et al. 1996) full grating scans of AB Aur 
were obtained in ISO revolutions 680 (at JD 2450717.747) and 
835 (JD 2450872.380), respectively. An SWS full grating scan 
of HD 163296 was made in revolution 329 (JD 2450367.398). 
Observing times were 3666 seconds for the SWS and 2741 sec- 
onds for the LWS observations. Data were reduced in a stan- 
dard fashion using calibration files corresponding to OLP ver- 
sion 7.0 (SWS) or 6.0 (LWS), after which they were corrected 
for remaining fringing and glitches. To increase the S/N in the 
final spectra, statistical outliers were removed and the detectors 
were aligned, after which the spectra were rebinned to a lower 
spectral resolution. The resulting spectra are shown in Fig. 1. 

A^-band (10.1 /im) images of HD 163296 were obtained 
on July 25, 1997 (at JD 2450654.608) using TIMMI on the 
ESO 3.6m telescope at La Silla. Total integration time was 
65 minutes. The pixel size was 07336, with a total field of 
view of 21' '5 x 21' '5. After a standard reduction procedure, 
the HD 163296 image was indistinguishable from that of the 
standard star r/ Sgr. After deconvolution of the HD 163296 im- 
age with that of r\ Sgr, the resulting stellar image stretched 
across 2 pixels. We conclude that the bulk of the 10 /im flux 
of HD 163296 comes from an area less than 077 (90 AU at 
122 pc) in diameter. 

3. Contents of spectra 

The infrared spectra of AB Aur and HD 163296 (Fig. 1) show 
big differences: whereas AB Aur shows the cool, strong con- 
tinuum expected for a Herbig star, for HD 163296 the contin- 
uum appears to be so weak that the entire SWS spectrum is 
dominated by solid-state emission features. The IRAS fluxes 
also plotted in Fig. 1 suggest that an underlying continuum 
in HD 163296 is present, but peaks longward of 100 /im and 
therefore is much cooler than the «40 K continuum in AB Aur. 



Table 1. Astrophysical parameters of programme stars 







Rpf 


1 1 \j 1UJZ7U 




a (2000) 


04 55 45.79 


(i) 


17 56 21.26 


(1) 


8 (2000) 


+30 33 05.5 


(i) 


-21 57 19.5 


(1) 


d [pc] 




(2) 


1221H 


(2) 


Sp. Type 


AOVe+sh 


(3) 


A3Ve 


(8) 


vn 


7.03-7.09 t 


(2) 


6.82-6.89 


(2) 


E(B - V) [ m ] 


0.16±0.02 


(4) 


0.04±0.02 


(4) 


Teff [K] 


r\i~f\f\-\-50Cl 

95001°°° 


(2) 


8700tjjjg 


(4) 


U [L ] 


47±12 


(2) 


26±5 


(4) 


R* [R©] 


2.5±0.5 


(4) 


2.2±0.2 


(4) 


M* [M ] 


2.4±0.2 


(2) 


2.0±0.2 


(4) 


log(Age) [yr] 


6.3±0.2 


(2) 


6.6±0.2 


(4) 


Ha 


P Cygni 


(5) 


Double-peaked 


(9) 


v sin i [km s _1 ] 


80±5 


(5) 


120±1 


(10) 


rot. period [ h ] 


32±3 


(6) 


35±5 


(11) 




76 


(7) 


58 


(7) 


P.A. [°] 


+79l 2 3 


(7) 


+1261 2 


(7) 



T However, older photographic measurements (Gaposchkin 1952) 
show values up to m pg = 8.4 for AB Aur. 

References: (1) Hipparcos Catalogue; (2) van den Ancker et al. 
(1998); (3) Bohm & Catala (1993); (4) This paper; (5) Bohm & Catala 
(1995); (6) Bohm et al. (1996); (7) Mannings & Sargent (1997); (8) 
Gray & Corbally (1998); (9) Pogodin (1994); (10) Halbedel (1996); 
(11) Catala et al. (1989). 

Both AB Aur and HD 163296 show a strong 9.7 /im amor- 
phous silicate feature in emission together with a broad emis- 
sion complex ranging from 14 to 38 /im. The emissivities for 
various dust components in the spectra are also included in 
Fig. 1 . In the HD 163296 spectrum the broad emission complex 
ranging from 14 to 38 /im is too broad and too intense to be 
solely attributed to the 19 /im feature due to the O-Si-O bend- 
ing mode. We tentatively attribute this feature to a blend of sil- 
icates and iron oxide. Lab spectra of FeO also show a strongly 
rising emissivity in the short-wavelength range of the SWS 
(Henning et al. 1995). When folded with a w 800 K black- 
body, this naturally produces a broad emission feature peaking 
around 3 microns, which is present in both stars. 

The large degree of redundancy in the SWS data makes it 
possible to assess the reality of weak spectral features which at 
first glance may appear to be lost in the noise. Each part of the 
spectrum was scanned twice by twelve detectors, so by check- 
ing whether a particular feature is seen in all detectors and in 
both scan directions, it is possible to disentangle real features 
from noise. The features identified in this way are listed in Ta- 
ble 2. AB Aur clearly shows the familiar 6.2 and 7.7 and pos- 
sibly also the 8.6 and 11.2 /im UIR bands usually attributed 
to emission by polycyclic aromatic hydrocarbons (PAHs), as 
well as a new UIR band at 15.9 /im. The 3.29 /im UIR band is 
absent. 

The HD 163296 spectrum shows a number of small emis- 
sion features at wavelengths corresponding to those of crys- 
talline olivines ((Mg x Fei_ x )2Si04). Remarkably, the amor- 
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Fig.l. ISO spectra of AB Aur (top) and HD 163296 (bot- 
tom) with the main features identified. The solid dots indicate 
ground-based and IRAS photometry (Berrilli et al. 1987, 1992; 
Hillenbrand et al. 1992; Weaver & Jones 1992). The dashed 
line in the top plot is the AB Aur spectrum after subtracting a 
spline fit to the continuum. Also shown (dashed lines at bot- 
tom) are the emissivities at a given temperature of the various 
dust components in the spectra (Jager et al. 1998; Dorschner et 
al. 1995; Henning et al. 1995; Bertie et al. 1969). The bars at 
the bottom of the plot show the IRAS 12, 25, 60 and 100 /jm 
passbands (FWHM). The inset in the AB Aur plot shows the 
UIR bands in the 5.5-9.0 /jm region after subtracting a pseudo- 
continuum. 



phous silicate feature has a higher band-strength relative to the 
continuum in our ISO data than in the ground-based 8-13 /xm 
spectrum of HD 163296 by Sitko et al. (1999). In contrast, the 
11.2 ^m shoulder, due to crystalline silicates, appears weaker 
in our spectrum, suggesting a significant time variability of 
both components. 

In addition to this, HD 163296 shows emission from the 
44 fim H 2 ice feature. The relative location of the IRAS 
60 [im measurement in comparison to the SWS spectrum sug- 
gests that the long-wavelength H 2 ice feature around 69 fim 
as well as the broad unidentified feature longward of 100 /zm, 
observed in HD 100546 and HD 142527 (Malfait et al. 1998, 
1999), might also be very prominent in HD 163296. PAHs are 
very weak or absent in HD 163296. 




log X [/Lim] 



Fig. 2. Spectral energy distributions of AB Aur (top) and HD 
163296 (bottom). Open squares and circles indicate observed 
and extinction corrected fluxes from literature (Wesselius et al. 
1982; The et al. 1985; Strom et al. 1972; Cohen & Schwartz 
1976; Berrilli et al. 1987, 1992; Hillenbrand et al. 1992; Weaver 
& Jones 1992; Mannings 1994; Gudel et al. 1989; Brown et al. 
1993; Skinner et al. 1993). The solid lines show Kurucz mod- 
els for the stellar photospheres of AB Aur and HD 163296. The 
dashed lines show linear fits to the submm data points for com- 
parison with the slope of the Rayleigh-Jeans tail of the Kurucz 
model. 

In addition to the solid-state features, both AB Aur and 
HD 163296 also contain a number of H I recombination lines at 
shorter wavelengths. All lines from the Bracket and Pfund in- 
cluded in the SWS wavelength range series are present, while 
the higher H I series are not detected, possibly due to the com- 
bined effect of a lower instrumental sensitivity and a higher 
background in this part of the spectrum. These recombination 
line data will be discussed in more detail in a forthcoming pa- 
per. 

The LWS spectrum of AB Aur is relatively smooth and fea- 
tureless. The only line that is clearly visible in the spectrum is 
the [Cll] line at 157.7 ^m. The strength of this line (8.6 x 
10~ 16 W m~ 2 ) is compatible with it originating in the back- 
ground rather than being circumstellar. As can be seen from 
Fig. 1, there is a «25% difference in the flux scales between the 
AB Aur SWS and LWS spectra in the overlapping region. Al- 
though within the formal errors of the absolute flux calibration 
for SWS and LWS, this discrepancy is larger than that found in 
other sources. The difference cannot be attributed to the dif- 
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Table 2. Contents of Spectra 



Feature 


FeO 


PAH 


PAH 


PAH 


PAH 


Si-O 


PAH 




PAH 


O-Si-O 


FeO 


A [fim] 


[3] 


[3.3] 


[6.2] 


[7.7] 


[8.6] 


[9.7] 


[11.2] 


[15.0] 


[15.9] 


[19] 


[23] 


AB Aur 


V 




V 


V 


V 


V 


v-- 


v-- 


V 


V 


V 


HD 163296 


V 


v-- 


v-- 






V 








V 


V 




Feature 


Oliv. 




Oliv. 


Oliv. 


Pyr.? 


Oliv. 


Oliv. 


Oliv. 




H 2 ice 




A [fim] 


[11.3] 


[13.4] 


[16.3] 


[17.8] 


[18.2] 


[23.5] 


[31.3] 


[33.5] 


[40.7] 


[43.8] 




AB Aur 














v-- 


v-- 


v-- 






HD 163296 


V 


V 


V 


V 


V 


V 


V 


V 


V 


V 





ferent aperture sizes (33" x 20" for SWS versus a circular 
80" FWHM for LWS), and confusion with extended emission, 
since then the LWS spectrum would have to have a higher flux 
level than the SWS spectrum. It is interesting to note that in 
the time interval between the SWS and LWS measurements, 
AB Aur did show an optical photometric event which could 
have also affected the infrared brightness (van den Ancker et 
al. 1999 and references therein), so the difference in flux be- 
tween SWS and LWS might in fact be due to real variability. 
In Fig. 1 we also plotted the IRAS fluxes of AB Aur. Although 
compatible with both spectra, the IRAS 60 fim flux agrees bet- 
ter with the SWS spectrum, so we rather arbitrarily choose to 
adopt the SWS flux calibration for the region around 45 [im. 

4. Spectral energy distributions 

Spectral Energy Distributions (SEDs) of AB Aur and 
HD 163296 were constructed from literature data as well as 
our new ISO spectra and newly obtained VLA photometry for 
HD 163296 and are shown in Fig. 2. All the submm fluxes in 
these SEDs refer to single-dish measurements. As can be seen 
from Fig. 2, the SED can be naturally decomposed in three 
parts: the optical wavelength range, where the total system flux 
is dominated by the stellar photosphere, the infrared to submm, 
where emission originates from the circumstellar dust disk, and 
the radio, where free-free emission from the stellar wind be- 
comes dominant. 

The difference in behaviour of the dust component in 
HD 163296 and AB Aur is striking: after a nearly flat en- 
ergy distribution in the infrared, the sub-mm and mm fluxes 
of AB Aur drop rapidly (\F\ oc A -4 3 ), indicative of the dust 
becoming optically thin at these wavelengths, whereas toward 
HD 163296 the slope of the sub-mm to mm fluxes (a A~ 29 ) is 
within errors equal to that of the Rayleigh-Jeans tail of a black 
body (oc A~ 3 ). The new radio points at 1.3, 3.6 and 6 cm for 
HD 163296 do not follow the simple power-law dependence 
expected if these were solely due to free-free radiation. This 
demonstrates that even at wavelengths as long as 1.3 cm, a sig- 
nificant fraction of the system flux is due to circumstellar dust. 
The 3.6 and 6 cm fluxes are probably dominated by free-free 
emission. 



The energy distribution of a circumstellar dust disk is gov- 
erned by its temperature profile, the density distribution and the 
dust properties (chemical composition and size distribution). 
Since the circumstellar disks of AB Aur and HD 163296 stars 
are expected to be passive (Waters & Waelkens 1998) and the 
properties of the central stars are nearly identical, the temper- 
ature profiles in the disks are expected to be similar as well. 
One possibility to explain the different sub-mm to mm slope for 
AB Aur and HD 163296 could be a much flatter density distri- 
bution for HD 163296. However, with a standard sub-mm dust 
emissivity (j3 = 2) the inferred dust mass for HD 163296 would 
become implausibly large. A better explanation may be that the 
dust properties of AB Aur and HD 163296 are different, a fact 
already concluded independently from the ISO spectra. To be 
able to radiate efficiently, dust particles must have a size similar 
to (or larger than) the wavelength, pointing to the existence of a 
population of mm- to cm-sized cold dust grains in the circum- 
stellar environment of HD 163296, whereas those in AB Aur 
must be micron-sized. The ISO spectrum of HD 163296 also 
contains warm (w 800 K) dust, suggesting a significant lack of 
emission from dust of intermediate temperatures. 

5. Discussion and conclusions 

We have shown that the main difference between the AB Aur 
and HD 163296 systems is that HD 163296 contains a pop- 
ulation of very large (mm to cm-sized), cold, partially crys- 
talline, dust grains, which is absent in AB Aur. AB Aur con- 
tains a population of small dust grains (PAHs), which is absent 
in HD 163296. In view of the fact that the stellar parameters are 
nearly identical except for stellar rotation, these differences are 
remarkable. This must mean that either the evolution of the dust 
composition in protoplanetary disks happens within the error in 
the age determination of both systems (2+ 2 Myr for AB Aur vs. 
4± 4 Myr for HD 163296; van den Ancker et al. 1998), or that 
the evolution of the dust is dominated by external factors. 

We have also shown that in AB Aur and HD 163296 iron 
oxide is a constituent of the circumstellar dust mixture and 
can also be responsible for the observed excess emission near 
3 /im. Since the near-infrared excess exhibited by our pro- 
gramme stars is by no means unusual for a Herbig Ae/Be star, 
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this means that the same applies for the entire group of Her- 
big stars. Therefore the results of models attributing this near 
infrared excess emission to very hot dust from an actively ac- 
creting disk (e.g. Hillenbrand et al. 1992) must be regarded 
with some caution. The case of HD 163296 demonstrates that 
infrared broad-band photometry can be completely dominated 
by emission from solid-state features, which must also be taken 
into account in any future modelling of the energy distributions 
of Herbig stars. 

The detection of PAHs in AB Aur shows that ground-based 
surveys (e.g. Brooke et al. 1993) have underestimated the frac- 
tion of Herbig stars containing PAHs. The case of HD 163296 
shows that the fraction of Herbig stars showing PAH emission 
will not go up to 100%, so models depending on the presence 
of very small dust grains to explain the observed near-infrared 
excess in Herbig Ae/Be stars (Natta et al. 1993; Natta & Kriigel 
1995) will not be successful in all cases. We believe iron oxide 
to be a more plausible explanation for this near-infrared excess. 

It is difficult to attribute the absence of the 3.29 fea- 
ture in AB Aur to an unusual temperature of the PAHs. Since 
the 6.2 and 7.7 /xm C-C stretches are strong, while the bands 
due to C-H bonds are weak or absent, a more promising possi- 
bility seems a very low hydrogen covering factor of the PAHs 
in AB Aur or the presence of a population of large (> 100 C 
atoms) PAHs (Schutte et al. 1993). If the new 15.9 urn UIR 
band in AB Aur is also caused by PAHs, this suggests that it is 
also due to a C-C bond. 

To gain more insight in the evolution of dust in protoplane- 
tary disks, it is useful to compare the spectra presented here to 
those of other Herbig Ae stars. Except for the absence of crys- 
talline material, the AB Aur spectrum and energy distribution 
are nearly identical to those of its older counterpart HD 100546 
(Malfait et al. 1998). In the case of AB Aur (2.5 M Q ), any 
possible crystallization of circumstellar dust must therefore 
occur at a stellar age older than 2 x 10 6 years. The differ- 
ences between HD 100546 and HD 163296 are larger: the crys- 
talline dust in HD 100546 is much more prominent than that in 
HD 163296 and the population of large cold dust grains seen in 
HD 163296 is absent in HD 100546. 

As these cases of AB Aur, HD 163296 and HD 100546 
demonstrate, the age of the central star and the degree of crys- 
tallization do not show a one to one correspondence, but the 
processes of grain growth and crystallization in protoplanetary 
disks are also not necessarily coupled. Studying a larger sam- 
ple of Herbig Ae/Be stars might shed more light on what causes 
this large observed diversity in dust properties in systems which 
appear very similar in other aspects. 
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